arXiv: 1507.08603v2 [cond-mat.str-el] 24 Nov 2015 


Expansion potentials for exact far-from-equilibrium spreading of particles and energy 


Romain Vasseur, Christoph Karrasch, and Joel E. Moore 
Department of Physics, University of California, Berkeley, CA 94720 and 
Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

(Dated: November 25, 2015) 

The rates at which energy and particle densities move to equalize arbitrarily large temperature and 
chemical potential differences in an isolated quantum system have an emergent thermodynamical 
description whenever energy or particle current commutes with the Hamiltonian. Concrete examples 
include the energy current in the ID spinless fermion model with nearest-neighbor interactions (XXZ 
spin chain), energy current in Lorentz-invariant theories or particle current in interacting Bose gases 
in arbitrary dimension. Even far from equilibrium, these rates are controlled by state functions, 
which we call “expansion potentials”, expressed as integrals of equilibrium Drude weights. This 
relation between nonequilibrium quantities and linear response implies non-equilibrium Maxwell 
relations for the Drude weights. We verify our results via DMRG calculations for the XXZ chain. 


The dynamics of how a system of interacting particles 
expands from an initial state with spatial variation of 
temperature, density, or both is one of the basic prob¬ 
lems in non-equilibrium statistical physics. The study of 
quantum effects on this process was reinvigorated by the 
experimental creation of ultracold atomic gases 00, in- 
eluding cases where the atoms are confined to one or two 
spatial dimensions. Originally the main quantity mea¬ 
sured was the momentum distribution 0®, but recent 
progress on the “quantum gas microscope” and related 
techniques has made it possible to image particle density 
with high resolution, e.g., on single sites of an optical 

lattice I5H3- 

Such imaging methods mean that important observ¬ 
ables to characterize expansion of an atomic gas in ei¬ 
ther free space or an optical lattice [8] are not the same 
as those for non-equilibrium processes in electronic trans¬ 
port. For electrons, the charge or energy current between 
two leads has been studied in hundreds of situations, in¬ 
cluding a few non-equilibrium results with interactions 
such as tunneling between Luttinger liquids [9] [10] , the 
interacting resonant level model Q3HH] and the single 
impurity Anderson model (for a recent review see [[15]). 
The point of the present work is to show that one nat¬ 
ural quantity of interest for atomic expansion measure¬ 
ments mm, namely the change in time of the first 
moment of particle or energy density, has a precise non¬ 
equilibrium thermodynamic description in a broad class 
of systems. For a continuum system with either Lorentz 
or Galilean invariance, this description reduces to stan¬ 
dard thermodynamic state functions, but we find that 
even lattice systems relevant to current experiments have 
a description in terms of an “expansion potential” that 
is distinct from conventional thermodynamic quantities. 

We use this description to compute the energy ex¬ 
pansion rate exactly in the anisotropic Heisenberg spin 
chain (XXZ model) and compare our results in detail 
against time-dependent density-matrix renormalization 
group (DMRG [22ti24] ) calculations using the finite tem¬ 
perature algorithm explained in [25] . The same for¬ 


malism is applicable to higher-dimensional systems with 
emergent Lorentz or Galilean invariance. Our predic¬ 
tions apply in particular to a one-dimensional Bose gas 
(Lieb-Liniger model [26j|27]> or its lattice regularization 
in terms of g-deformed bosons [28]). expanding into vac¬ 
uum, a problem that has attracted a lot of attention re¬ 
cently [81 129H40] . Our results show that, at least for some 
quantities, exact results can be obtained for far-from- 
equilibrium expansion even in lattice models at arbitrary 
coupling strength. 

At t = 0, prepare two semi-infinite regions x < 0 and 
x > 0 at equilibrium with chemical potentials and tem¬ 
peratures and (iir,Tr) (Fig [lji) . (The initial 

state on the boundary between the two leads, or a possi¬ 
ble finite extent of the boundary region, will not matter 
for the quantities of interest here after some initial tran- 
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FIG. 1. (a) Nonequilibrium expansion setup considered in 

this letter, (b) Energy point current jE(x,t) in the XXZ spin 
chain (see Fig [5] for parameters), (c) The variation of the 
expansion potential Q does not depend on the path in the 
(/ 6 , fi) space. This implies nonequilibrium Maxwell relations 
(see text). 
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FIG. 2. Spatially integrated charge and energy currents in the 
XXZ model with open boundary conditions as a function of 
time for one choice of parameters Tl, /xl, Tr, / lr at A = 0.5 
under two protocols (A) and (B) that differ only in the way 
the central bond is dealt with in the initial state [JJ. The 
currents J — f^ A jdx are locally integrated around the cut 
site. The chemical potentials pl,r prepare the state but are 
not included in the real-time evolution (i.e., they are chemi¬ 
cal potentials rather than electric potentials). Inset: energy 
density profile as a function of time. The spatially integrated 
energy current is equal to dMl h /dt where Ml h is the first 
moment of energy density. 


sient.) We write one-dimensional equations for simplicity 
but the concept is general. We quantify the expansion 
for t > 0 by the time dependence of the first moment of 
particle density, or similarly for energy, 

Mi(t) = f n(x,t)xdx , (1) 

J -A 

with A a large observation scale: A vt with v a typical 
velocity. From now on we suppress the arguments of n 
and Mi. The continuity equation relates density and 
current d t n + d x j = 0. Now = — f^ A xd x j dx = J, 

with J = f^ A j dx where in the integration by parts we 
have assumed j(x) « 0 at x = ±A (Fig[^>). 

The key ingredient for the existence of an expansion 
potential is the conservation of integrated current: 


j dx , H 


= 0 , 


( 2 ) 


which is true for many problems of interest with periodic 
boundary conditions. Note that this is a stronger state¬ 
ment than what is sometimes meant by a “conserved cur¬ 
rent” , which is anything related to a conserved charge by 
a continuity equation. A simple example with such a con¬ 
servation law is a Bose gas in d spatial dimensions with 
say, ^-function interactions H = J d d x + 
with [\I/f (x), 4/(x/)] = S(x — y ), where the to¬ 


tal particle current Jq = —i J — is con¬ 

served. More generally, a system with one species of par¬ 
ticles moving in the continuum in any spatial dimension 
will satisfy © for particle current if particle current is 
proportional to total momentum and momentum is con¬ 
served by the interactions. A less trivial example of (J2| 
is energy current in the spinless fermion model or XXZ 
spin chain (we will use the former representation): the 
energy current operator Je = i hj+i] commutes 

with the XXZ Hamiltonian H xxz = JT hi with 


, J 
hi = ~2 


c \+i c 


J A 




(3) 

with rii = cjci, implying purely ballistic energy trans¬ 
port in 02 ]. An example of a current not conserved 
in this sense is charge current in the XXZ model; while 
there is a degree of ballistic transport in this model in the 
gapless regime even at nonzero temperature [43H45] . the 
commutator in © is nonzero. Steady-state energy cur¬ 
rents between reservoirs have been actively studied [Mb 
[52] but exact results have been difficult to obtain except 
in the low-temperature conformal limit or for noninter¬ 
acting systems. 

Expansion potentials. The global current conserva¬ 
tion law © implies that the current density should itself 
satisfy a continuity equation for some “current of cur¬ 
rent” P, 


d t j + d x P = 0, (4) 

and we will see in the following that the operator P is 
related to pressure for systems with emergent Galilean or 
Lorentz invariance. Now spatially integrate this second 
continuity equation © over the region [—A, A] centered 
on the boundary between our two large reservoirs L and 
R. Then 

= -P\t a = a s = g L - s R , ( 5 ) 

where we have introduced the expansion potential 
G(fj,,T) = (P)h,t fc> r the thermodynamic expectation of 
the operator P. 

This is a strong constraint on the integrated current 
J = J_ A jdx. We perform DMRG calculations on a XXZ 
spin chain with open boundary conditions that effectively 
describes a region of an infinite system. Within that re¬ 
gion (shown in Fig.JTJ)), the total energy current is clearly 
not conserved [41 and grows linearly with time (Fig. [2| 
for times short enough that the reservoirs are effectively 
infinite, so their initial values can be used in the bound¬ 
ary evaluation on the right-hand side of ©• If the cur¬ 
rent has both diffusive and ballistic components (like the 
charge current in the XXZ chain), diffusive contributions 
die out after a transient and the spatially integrated cur¬ 
rent also grows linearly. However, the situation becomes 
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FIG. 3. Comparison between charge Drude weight (long-time 
asymptote of (J7 q(£)J7q( 0))/LT and rate of particle spread¬ 
ing d 2 Mi/dt 2 = (J7q (£))*/£ for a small chemical potential 
difference (A fi ~ 10 -3 <C J = 1) in the XXZ chain (pro¬ 
tocol A). Inset: similar relation for energy transport be¬ 
tween thermal Drude weight and rate of energy spreading 
d 2 M{ h /dt 2 = (sJE(t)) t /t for a small temperature difference at 
half-filling (see eq. dm)). 


especially simple for a current satisfying ©, the key be¬ 
ing that the right-hand side of (J5| contains only the op¬ 
erator P evaluated at equilibrium , since deep within the 
reservoirs the system remains arbitrarily close to equilib¬ 
rium in this intermediate time regime. This result relies 
only on © and does not depend on whether the system 
is gapped or gapless for instance. 

Linear response. One more relation is all that is 
needed to compute the expansion potential in some im¬ 
portant cases. This is because eq. © implies that linear- 
response is enough to predict non-equilibrium, since lin¬ 
ear response gives the derivative of Q , and knowing its 
derivative determines the function up to an arbitrary 
additive constant. Focusing for the moment on en¬ 
ergy current and a purely thermal gradient, linear re¬ 
sponse then predicts j E = — cr E NT with the thermal 
conductivity characterized by a thermal Drude weight 
c t E (u ) = 7rD th (T)S(uj) with Ah = P 2 (J e )/L, where 
L is the size of the system. The spatially integrated 
current between the two reservoirs R and L then reads 
f- A j E dx = ttAT5(uj = O)Ah(T) where the time t 
can be thought of as an infrared cutoff that regularizes 
S(cu = 0) ~ f l _ t H =t/ 7r. We thus find 

(R Tlfth i 

^^ = -(JE)t=Ah(T)x(AT), (6) 

with J E = f^ A j E dx and (... )t refers to the nonequi¬ 
librium expectation value after time t. For the charge 
current at constant temperature Tr = T E = T, we sim¬ 
ilarly find d J^ 1 = (— D c A/a with the charge 


Drude weight D c = /3(Jq)/L (if [H, Jq\ = 0), for a small 
chemical potential gradient A (i. These results are easily 
extended to the case where both temperature and chem¬ 
ical potential gradients are present (see below). We also 
note that these linear response results remain valid even 
if the currents are not fully conserved and contain diffu¬ 
sive parts, like the charge current in the XXZ spin chain, 
which provides a direct way to measure Drude weights 
via imaging in cold atom experiments (see also [53]). 
We checked this relation between charge ( resp . thermal) 
Drude weight and linear-response rate of spreading of 
charge (resp. energy) in the XXZ chain (see Fig. [3]) - 
similar relations also exist for diffusive systems [54] . 

Nonequilibrium expansion potentials. The thermody¬ 
namic description eq. © together with the linear re¬ 
sponse prediction implies that the spreading of particles 
and energy far from equilibrium are fully characterized 
by the equilibrium Drude weights. As an example, let us 
consider the rate of energy spread in the XXZ spin chain 
between two reservoirs at different temperatures Tr and 
T e and fi = 0. Then even far from equilibrium 

dM th C Tl 

In other words, the nonequilibrium rate of energy spread 
is given by the variation A r^lGe = Ge(T l ) ~Ge(Tr) of 
a state function Ge(T ) with OtGe = Ah(T). This can be 
checked numerically by comparing the rate of expansion 
to the thermal Drude weight of the XXZ model computed 
by Klumper and Sakai ES] (see Fig. [I). 

This is easily generalized to the case of reservoirs R 
and L with both different temperatures (Tr and T E ) and 
chemical potentials (/ lr and /ul) . If the energy current is 
conserved, eq. © implies that the far-from-equilibrium 
rate of energy spread is given by the variation of an ex¬ 
pansion potential Ge(p,P = T -1 ) 


d 2 M th C 

i,<2 Ar^lGe = / dGE , ( 8 ) 

dt Jr^l 

where the differential dGE is exact so that the integral 
does not depend on the chosen path. The state function 
Ge is then fully determined by the equilibrium Drude 
weights associated with the conservation of the energy 
current. Linear response theory [56] then yields 


d g E = pJMA du 



dp. 


( 9 ) 


Even if Jq is not conserved, the Drude thermopower is 
a thermodynamic quantity determined by ( JqJe) pro¬ 
vided that {Jee H] =0. If the particle current is con¬ 
served, we find similarly that the integrated nonequi¬ 
librium particle current between two reservoirs (priPr) 
and (pll,Pl) is given by the variation of another state 

function = A r^lGq with dQ Q = /3^p-dn - 

( {JqJe)_ _ ^ Vspjdp. 
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FIG. 4. The thermodynamic description § implies that the 
spatially integrated current Ji ^.2 between two reservoirs 1 
and 2 should be equal to Ji _>3 + for any intermediate 

reservoir 3. We verified this “cyclic invariance” of the spa¬ 
tially integrated energy current (jE(t))t/t and point current 
(j e( x — 0 jt))t in the XXZ chain (protocol A). While cyclicity 
for the point current may be only approximate, it is exact for 
the integrated current. Insets: numerical check of eq. 0 and 
of the nonequilibrium Maxwell relation (10) that follows from 
conservation of energy current. 


Nonequilibrium Maxwell relations. We saw above 
that when either the energy or particle current is fully 
conserved, then even far from equilibrium the expan¬ 
sion dynamics of energy or particle densities are charac¬ 
terized by state functions that are entirely determined 
by equilibrium Drude weights. An interesting corol¬ 
lary of the path-independence of these state functions 
(Fig §) are nonequilibrium Maxwell relations for the 
Drude weights. For example, if the energy current is 
conserved, O^O^Qe = dp O^Ge yields 

(fidp ~ JqJe) + Q^Je) = 0? (10) 

which can also be rewritten as (A HJqJe) = (ANJ^) 
with AH = H - ( H) and AN = N - (N). This equal¬ 
ity was known in the context of the XXZ chain m and 
was actually used to compute the Drude thermopower 
analytically [58], but our approach provides a very trans¬ 
parent derivation of why such a relation has to hold (see 
Fig. 0 for a numerical check). If the charge current is 
conserved, then the associated nonequilibrium Maxwell 
relation reads (ftdp - pbd tJL )(J q) + Q^JeJq) = 0, which 
can also be rewritten as (A HJq) = (A N JqJe)- 

Examples in d > 1 dimensions. Even though most 
of the arguments discussed above focused on one dimen¬ 
sion for simplicity, the general concepts apply in higher 
dimension as well. For a system with emergent Lorentz 
symmetry (z = 1 critical points for instance), the symme¬ 
try of the stress-energy tensor means that the energy cur¬ 
rent Ten with i = 1,..., d is also the (conserved) momen¬ 
tum density T^q. The energy expansion potential then 


reads Ge{P) = ~ F dp f d d x^Ei( T oi(x)T Oi (0)), which 
can be related to pressure ISBE2I. In a non-relativistic 
system with a single species of particles and current pro¬ 
portional to (conserved) momentum, there is a particle 
expansion potential; the interacting Bose gas is one such 
example. The particle Drude weight D c is then entirely 
determined by the sum rule f ^a(ou) = D c = — with 
n the density and m the mass [56] . This immediately 
implies that the expansion potential is simply related to 
pressure Qq = — = pp with Q the thermodynamic 

grand potential and V the volume - this is a consequence 
of Galilean invariance m- The Drude thermopower is 
then given by {JqJe)/V = ^(u + P) with u the inter¬ 
nal energy density. These quantities can be computed 
explicitly for the Lieb-Liniger gas in one dimension as a 
function of T and pi (or particle density) [60] . This and 
other simple cases where the expansion potentials can 
be computed explicitly, such as non-interacting systems 
and Luttinger liquids, are given in Supplemental Mate¬ 
rial [4T. 

Nature of the steady-state. Interestingly, the varia¬ 
tion of expansion potential A Q provides a lower bound 
for the point current j(x) [50]. However, the more general 
relation between spatially integrated and point currents 
remains mysterious. We find numerically that both the 
energy density ue{xJ) and the energy current j^(x,t) 
in the XXZ spin chain at half-filling become functions 
of x/t at large enough times, with nontrivial limiting 
shapes mt In the low-temperature limit described by 
conformal field theory 03 S3, we expect a uniform 
steady-state local current Je(x) = ^§y = — Tj|) 

over a region of size 2 vt with v the spinon velocity. How¬ 
ever, we find that the rescaled functions jE(x/t), nE{x/t ) 
even at moderate temperatures are very far from that pic¬ 
ture: in general, there is no nonzero range of the reduced 
variable x/t for which jT(x/t) is constant, indicating 
that the steady-state region spreads sub-ballistically, and 
there are no transient “shock-waves” like those expected 
in the presence of Lorentz invariance 152 EH separating 
the uniform steady-state region from the reservoirs. It 
is an interesting problem for future work to determine 
more properties of the limiting function jE(x/t), pos¬ 
sibly by adapting the recently developed hydrodynamic 
approaches for relativistic systems [52] [61] to incorporate 
the additional conserved quantities of integrable lattice 
spin chains. 

Discussion. In closing, we emphasize that the expan¬ 
sion potentials generalize familiar concepts in the pres¬ 
ence of either Galilean or Lorentz invariance to consid¬ 
erably more complex physical situations. Lattice models 
for which a current is conserved in the sense of § in¬ 
clude the XYZ spin chain, the g-Bose gas [28], and the 
supersymmetric point of the t-J model [62] . For systems 
where the conservation law does not strictly hold, such as 
the Bose-Hubbard model at small occupancy where rare 
double occupancies spoil the mapping to the XXZ model, 
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Joule heating and other strongly non-equilibrium physics 
could be computed using perturbation theory from the 
expansion-potential case. It would be interesting to con¬ 
nect the expansion potential to other nonequilibrium ef¬ 
fects, such as “quantum quenches” of a coupling 63] . 
which can reveal topological phases [64, 65 . For lat¬ 
tice models with conserved energy current but without 
full integrability, the expansion potential still exists and 
could be computed numerically at equilibrium , while it 
would serve as a useful constraint on predictions about 
far-from-equilbrium energy flow m- 
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NUMERICAL METHODS AND BOUNDARY CONDITIONS 

We run our DMRG calculations using the following numerical parameters. The system size is typically chosen as 
L = 200, and for a few cases we check our results against L = 100 and L = 300 to ensure that we are effectively in the 
thermodynamic limit. The imaginary time evolution from T = oo down to the desired temperature is carried out using 
a second order Trotter decomposition with a step size of Sr = 0.01, and the discarded weight during each individual 
sub-step is at most 10 -12 . The real time evolution is performed using a fourth order Trotter decomposition with a 
step size of St = 0.2. The calculations are performed for various values of the discarded weight, which eventually is 
chosen such that physical quantities have a relative accuracy of approximately 1 percent. Typical values are 10 _7,_9 
for the linear-reponse correlators as well as 10 _9,,,_n for the non-equilibrium currents at small bias in Fig. 2 or 10 -6 
for the currents in presence of a large bias in Fig. 3. A simple test of the method is that the free case is reproduced 
exactly and the results do not change on increasing the size of the matrix product states used in the calculation. 

We use numerically two protocols A and B to prepare the XXZ spin chain out of equilibrium. For protocol A, the 
system is prepared with the density matrix e~dL(H L -ULN L )-f3 R (H R -^ R N R ) ^ where Hl and Np are the Hamiltonian and 
number operator in the left system including the bond to the right system. The chemical potential at the interface 
(i.e., at the first site of the right system) is chosen as (ul + Ur)/ 2. For protocol B, this central bond is initially cut, 
and the density matrix reads e~0L{H L -^ L N L ) ^ e (3r(h r n R N R ). The chemical potentials Ul,r prepare the state but 
are not included in the real-time evolution (i.e., they are chemical potentials rather than electric potentials). We find 
that both protocols give the same long-time limit for the spatially integrated currents Je/ t and J7 q/£. 

We now comment on the boundary conditions used in our DMRG calculation, and explain how an exact conservation 
law is consistent with the change in the integrated current observed in our numerics. As discussed in the main text, 
the conservation law [H , Je\ = 0 for the XXZ spin chain is valid only for a system with periodic boundary conditions, 
with Je = § jEdx. For such a periodic chain with two halves A and B prepared at different temperatures Ta and 
Tg, our main results apply to the partially integrated current J^ A jEdx = t x (Qa — Gb) centered over one of the two 
junction regions between A and B. Note that this does not contradict the conservation of the total spatially integrated 
current, as there is a net partially integrated current moving in the opposite direction at the opposite junction, so that 
the total current f jEdx = t x (Qa ~ Gb) + t x (Q B — Qa) = 0 as required by the conservation law (Fig. 1). Our DMRG 
method is applied to a XXZ spin chain with open boundary conditions that effectively describes a region of an infinite 
system, and clearly within that region the total current is not conserved, but there are non-trivial consequences of 
the local current conservation OJe + 9 x Pe = 0 as described in the main text. 


BALLISTIC TRANSPORT 

In this appendix, we provide a brief review of ballistic heat transport for the sake of completeness. Perhaps the 
most convincing experiment showing an extreme violation of Fourier’s law for heat transport is the observation of 
the “thermal conductance quantum” [1]. We will review this famous ID example of how heat in certain systems is 
characterized by ballistic transport and a quantized conductance rather than a finite conductivity a la Fourier. The 
XXZ model analyzed in the main text, which appears frequently in atomic physics, is an example of an interacting 
model that, because of its exact integrability, shows ballistic properties in its heat conduction, but it seems worthwhile 
to first discuss the free case before moving to the non-trivial interacting case. 

What is the thermal conductance of a free gas of one-dimensional particles, e.g., bosons such as the phonon? One 
might think it is infinite as they propagate without scattering, but as for electrons it pays to be careful. The charge 
conductance of a ID gas of electrons at low temperature was shown to be G = e 2 /h “per channel”, i.e., per Fermi 
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FIG. 1: (a) Energy transport of free particles: the left and right halves independently radiate at different temperatures, and the 
difference in radiated energy determines the thermal current, (b) Segment of an infinite ID system with left and right halves 
initially prepared at different temperatures. Energy density evolves as shown in the region of the boundary, which means that 
there is a net rightward current of energy. This obviously does not violate conservation of spatial integral of energy current-the 
same local increase would hold for a segment of a large periodic system until such times as the two junction regions between 
high and low temperatures notice each other. In an infinite system, that time diverges. 


surface point, rather than being infinite. What about the energy current if the two leads are at different temperatures 
instead of different voltages? 

We can use the Landauer approach to compute this: for the case of perfect transmission, this amounts to computing 
the right-moving energy current from a lead at temperature T + dt and subtracting the left-moving energy current 
from a lead at T. Assume one-dimensional free bosons as in the Schwab et al. experiment mentioned above. Using k 
for momentum, we have that the total energy current (units of energy per time) is 


Here v & 
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Here cj max is the highest phonon frequency. If we assume that the temperature is small compared to this, so that 
x = huj/ksT runs from 0 to infinity, then we obtain (note that we need to multiply by ( ksT/h ) 3 ) 


Je 



dx 


x 2 e x 

(e x - l) 2 ' 


( 3 ) 


The dimensionless integral gives 7 t 2 /3, so 


= Je = EkET 
dt 3 h 


( 4 ) 


An interesting fact about the thermal conductance Go is that it is the same for bosons or fermions (or indeed anyons), 
unlike charge transport. The Schwab et al. experiment observed one thermal conductance quantum Go for each 
low-temperature phonon mode. 

Ballistic heat transport consistent with (4), with a factor of central charge, was recovered in the low-temperature 
limit of any ID CFT by Sotiriadis and Cardy [2] and by Bernard and Doyon [3]. This result can be viewed as a 
form of the Stefan-Boltzmann law (see Fig. la): when particles are moving freely, the left half-line and right half-line 
both radiate particles that move directly through each other, so the thermal current is related to the difference in the 
Stefan-Boltzmann law of total radiated power between the left and right halves. 
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FIG. 2: Expansion potentials Gq(/jl,/3) (a) and Ge(h,P) (b) characterizing the spatially integrated nonequilibrium spin and 
energy currents in the XX chain. 


EXPANSION POTENTIALS FOR NON-INTERACTING FERMIONS 


In this appendix we compute the energy and spin/charge expansion potentials Qq(P,h) and Ge(P^) in the XX 
chain (XXZ at A = 0) 

ffxx = j X (SfSf +1 + W+i), (5) 

i 


which can be mapped onto non-interacting fermions with dispersion e/~ = — cos k — here and in the following, we 
will set J = 1 . Both the charge current Jq = ^ k VkTik and energy current Je = Yh k v k^k^k commute with the 
Hamiltonian Hxx = ^k e k n k, where n k = c\c k is the occupation number of the level k and Vk = d^k = sinfc. 
Therefore, as explained in the main text, the nonequilibrium spreading of particle and energy are fully characterized 
by the equilibrium Drude weights (J7j), (Jq) and ( JeJq )• Those can readily be computed from the generating 
function 

g(P,n ,\ Q ,\ E ) = -y\nTr e-^ H -^ N)+x ^ J ^ +XEjE . ( 6 ) 

±J 


It is then straightforward to verify that they satisfy the nonequilibrium Maxwell relations (ft dp — /x3 M )( JqJe) + 
®h(Je) = 0 and (/ 3dp — fid^) (Jq) -*r O^JeJq) = 0 as they should. We also remark that since both charge and energy 
currents are conserved, the charge and energy grand expansion potentials are related through O^Qe — (/iS^ — /3dp)GQ- 

The expansion potentials are then obtained by integrating the differentials (IQe = P dfi — — fi ^ d(3 

and dg Q = d p. 


and 


Gq(P^) 


Qe(P^) 


r dk vi 

L 27r 1 + ’ 

r dk v k e k 

J_ n 27r 1 + ' 


( 7 ) 

( 8 ) 


Focusing on energy transport, Ge(P^) = — J + ^T 2 + ... at half-filling and low temperature. We remark that the 
T 2 term is entirely fixed by conformal field theory: it coincides with 2 VEg(T) where ve = 1 is the Fermi velocity 
and g(T) = jt>T 2 is the function characterizing the energy point current Je — q(Tl) — q(Tr) in the steady-state 
region [3] (and we have used the fact that the central charge c — 1 for the XX chain). This is consistent with a 
steady-state region expanding as [—vpt] with uniform point-current Je = q(Tl) — #(Tr), leading to the spatially 
integrated current f jEdx = 2 VFt(g(TL) ~ q(Tr)) = t^GiTj) — G(Tr)). We emphasize however that this relation 
between point-current and expansion potential holds only in the low-energy regime. Whereas we have argued that 
the far-from-equilibrium spatially integrated current can be computed (even at high temperature) in terms of an 
expansion potential fully characterized by equilibrium Drude weights, it remains unclear whether a similar relation 
also holds for the point current. 
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FIG. 3: Grand expansion potential Gq(HiP) = ^ characterizing the non-equilibrium rate of particle density spread in the 
ID Bose gas with (5-function interactions with interaction strength c — 1. Inset (a): influence of the interaction strength on 
Qq(HtP) and comparison with the c = 0 (non-interacting Bose gas) and c = oo (effectively free-fermions) limits. Inset (b): 
Thermopower Drude weight determined from the nonequilibrium Maxwell relation. 


EXPANSION POTENTIALS FOR A LUTTINGER LIQUID 

In this appendix, we derive the form of the particle and energy expansion potentials for a Luttinger liquid. The 
low-energy degrees of freedom of the XXZ chain in the massless regime can be described by the Luttinger Hamiltonian 

H=^Jdx ( K(d x 9 ) 2 + K-\d x< j>f) , (9) 

where the velocity v and the Luttinger parameter K depend in a non-universal way on microscopic parameters. The 
chemical potential /i is then introduced by coupling to the particle number N = f dxd x (j)/^K - note that we consider 
v and K to be constants in the field theory description, even though they depend on the chemical potential/Fermi 
velocity of the underlying microscopic model. The particle and energy currents read Jq = vK f dxd x 6 / y/ir and 
Je = v 2 f dxd x 0d x (j) ; they both commute with the low-energy Hamiltonian H. The corresponding Drude weights are 
easily computed [4, 5] - note that in particular, ((Je ~ hJq)Jq) = 0. This yields 

dg Q = —d^, (io) 

IT 

IXVK 7TV , . 

dg E = ~E dli ~ 3p3 d P- ( n ) 

Within this low-energy limit, the expansion potentials are simply related to the value of the steady-state local (point) 
currents jq = ^(/^l — Hr) and Je = j — T^) + ^(Hl ~~ Hr) by a factor 2v (the size of the steady-state region 
being ~ 2 vt). 


PARTICLE DENSITY EXPANSION IN THE INTERACTING BOSE GAS 


In this appendix, we discuss the particle expansion dynamics in the interacting Bose gas with (5-function interactions. 
We start with the one-dimensional (integrable) Lieb-Liniger model 

ff = -X^+ 2 c X^-^)- ( 12 ) 

i 1 i<j 

The Lieb-Liniger model was realized experimentally in various circumstances [6-9]. Since the particle current is 
conserved, the far-from-equilibrium expansion dynamics is entirely determined by the state function Qq. The 
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T l =°°, Tr=1/2 


T l =1 , T r =1/8 



(a) 

T L =1/2, T r =1/4 



x/t 

(c) 

FIG. 4: Energy point current jE(x,t) as a function of x/t i 
parameters Tl and Tr. 



(b) 

T L =1/4, T r =1/8 



x/t 

(d) 

the XXZ spin chain at half-filling for various values of the 


Drude weight ( Jq) /L can be obtained from the generating function — /3g(/3, /x, A) = In (Tr e -P( H ~v N )+ x JQ^ as 
(Jq)/L = d\(—f3g) | . Since the Lieb-Liniger model is integrable, the function g(/?,/x, A) can be computed using 

the Thermodynamic Bethe Ansatz (TBA) [10, 11] as g((3, /x, A) = — T f ln(l + e _e ^^/ T ) with e(k) the pseudo-energy 
dressed by the interactions that can be determined from the non-linear integral equation 

e(k)/T = (3(k 2 - /x) - 2Xk - 1 J g2 + ° _ In (l + e~< k ^ T ) dq. (13) 

The nonequilibrium Maxwell relation (/3dp — /x3 m )(J7q)+S m (J7£,7q) = 0 then gives the Drude thermopower ( JeJq)/L , 
which is an interesting quantity since thermoelectric properties can be measured in atomic gases [12], and the nonequi¬ 
librium Maxwell relation is crucial to compute it. Here we are working at fixed chemical potential rather than number 
density, but note that as the expansion process takes place in a closed quantum system and the initial reservoirs are 
thermodynamically large, the expansion rate for a density difference follows from working at the chemical potentials 
corresponding to the initial densities. 

As mentioned in the main text, the particle Drude weight can also be obtained much more easily from the sum 
rule f ^f-cr(cj) = /3 ® ^ with n the density and m = \ the mass. This means that expansion potential is simply 

related to pressure Qq = — ^ with D the thermodynamic grand potential. The Drude thermopower is then 

given by (JqJe) /L = ^(u + P) with u the internal energy density. We plot these quantities as a function of /3 and g 
for the Lieb-Liniger gas with fixed interaction strength c = 1 by solving numerically the TBA equations (see Fig. 3). 
We emphasize that the relation Qq = P/m between the expansion potential and pressure is general and does not 
rely on integrability as it follows from Galilean invariance, and also holds in higher dimensions. The only step where 
the integrability of the Lieb-Liniger model was crucial is to compute explicitly P as a function of temperature and 
chemical potential or particle density. 
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Even though the particle expansion potential becomes especially simple in the continuum limit for systems with 
emergent Galilean invariance, we emphasize that its existence is more general. In particular, the expansion potential 
Qq can also be computed exactly for the so-called q -Bose model [13], one of the (integrable) lattice regularizations of 
the Lieb-Liniger model. The q-Bose model has an exact conserved particle current not related to momentum, and the 
corresponding Drude weight can be computed using TBA equations [11]. The Drude thermopower and the expansion 
potential Qq then follows from integrating the differential cIQq. 


NATURE OF THE NONEQUILIBRIUM STEADY-STATE 

To investigate further the nature of the steady state and transient regions, we plot the point energy current t) 
against the reduced variable u = x/t in the XXZ chain at A = 0.5 for various temperatures T^, Tr and fi = 0. We 
find that the current converges to a limit shape Je = f{x/t\ with f duf(u) = A Q. At low temperature, conformal 
field theory results imply the existence of a ballistic steady-state region with uniform current [3, 14], spreading at 
the sound velocity v, sharply separated by “shock waves” [15, 16] from the rest of the system outside the light 
cone. Our numerical results (Fig. 4) suggest that even at moderate temperatures, the irrelevant perturbations to 
the Luttinger liquid theory breaking Lorentz invariance modify that picture dramatically, reducing the steady-state 
region to u = x/t = 0 - indicating a steady-state spreading sub-ballistically. We also observe broad transient regions 
that propagate linearly with time. In the non-interacting case (XX spin chain), the limit shape can be expressed as 
JE = f(x/t) = g(x/t,T L ) - g(x/t,T R ), with 

g{u = x/t,T) = / 1 EEJt (% + v k) ~ °( u ~ v k )) i ( 14 ) 

with 6 the Heaviside step function, e& = — cos k and Vk = dk e k- The interacting case is much more involved, and 
it would be very interesting to investigate whether the recently developed hydrodynamic approaches for relativistic 
systems [15, 16] can be modified to describe such nonequilibrium steady-states in integrable spin chains. 
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